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Abstract 

A study of the inclusive charged hadron production in two-photon collisions is 
described. The data were collected with the DELPHI detector at LEP II. Re- 
sults on the inclusive single-particle px distribution and the differential charged 
hadrons da/dpx cross-section are presented and compared to the predictions of 
perturbative NLO QCD calculations and to published results. 
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1 Introduction 

The inclusive production of hadrons in 7*7* interactions can be used to study the 
structure of two-photon collisions [1]. These photons are radiated by beam electrons 
which scatter at very small angles and most of them are not detected. The untagged 
photons are quasi- real with a mass Q 2 ~ 0. At LEP II these collisions are the main 
source of hadron production, providing a good opportunity for such an investigation and 
thus to check the predictions of leading and next-to-leading order (NLO) perturbative 
QCD computations. 

The L3 and OPAL collaborations have published results of their analyses of the inclu- 
sive production of charged hadrons in two-photon collisions [2,3]. While L3 observes a 
pion production cross-section largely exceeding the NLO QCD predictions at high trans- 
verse momenta (5 GeV/c < p T < 17 GeV/c), OPAL finds a good agreement with them, 
in the Pt < 10 GeV/c range of its analysis. 

In this paper we present the DELPHI study of the inclusive production of charged 
hadrons in collisions of quasi-real photons. Section 2 describes the selection criteria 
for the event sample collected for this study. The inclusive single-particle transverse 
momentum spectrum and the measurement of the differential charged hadrons cross- 
section are presented in Section 3. They are compared to theoretical QCD predictions 
and published results in Section 4. 

2 Experimental procedure 

The analysis presented here is based on the data taken with the DELPHI detector [4,5] 
in 1996-2000, covering a range of centre-of-mass energies from 161 GeV to 209 GeV, with 
a luminosity-weighted average centre-of-mass energy: 195.5 GeV. The selected data set 
corresponds to the period when the Time Projection Chamber (TPC), the main tracking 
device of DELPHI, was fully operational thus ensuring good particle reconstruction. The 
corresponding integrated luminosity used in this analysis is 617 pb _1 . 

The charged particles were measured in the tracking system of DELPHI, which consists 
of the micro Vertex Detector (VD), the Inner Detector (ID), the TPC, the Outer Detector 
(OD) in the barrel, and the Forward Chambers FCA and FCB in the endcaps of DELPHI, 
all embedded in a homogeneous 1.2 T magnetic field. The following selection criteria are 
applied to charged particles : 

• transverse momentum pt > 150 MeV/c; 

• impact parameter of a trajectory transverse to the beam axis A xy < 0.4 cm; 

• impact parameter of a trajectory along the beam axis A 2 < 2 cm; 

• polar angle of a track with respect to the e~ beam 10° < 9 < 170°; 

• track length / > 30 cm; 

• relative error of its momentum Ap/p < 100%. 

The measurement of neutral particles is made using the calorimeter information pro- 
vided by the electromagnetic calorimeters, the High Density Projection Chamber (HPC) 
in the barrel and Forward Electromagnetic Calorimeter (FEMC) in the forward (back- 
ward) regions and by the hadronic calorimeter (HAC). Events with photons tagged by 
the DELPHI luminometer (STIC), i.e. with high Q 2 values, have been rejected. The 
calorimeter clusters, which are not associated to charged particle tracks, are combined to 
form the signals from the neutral particles (7, ir°, K° L , n). The following thresholds are 
set on the measured energy: 0.5 GeV for showers in the electromagnetic calorimeters and 
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2 GeV for showers in the hadronic calorimeter. Furthermore the polar angle of neutral 
tracks was required to be in the range 10° < 9 < 170°. 

To extract the hadronic events from the collisions of quasi-real photons the following 
cuts are applied: 

• energy deposited in the DELPHI luminometer (STIC: 2.5° < 9 S tic < 9°) 
Estic < 30 GeV; 

• number of charged-particle tracks N c h > 4; 

• visible invariant mass, calculated from the four-momentum vectors of the mea- 
sured charged and neutral particles, assuming the pion mass for charged particles, 
5 GeV/c 2 < W vis < 35 GeV/c 2 . 

The first condition eliminates the so-called single and double-tagged 7*7* events. The 
condition on the charged track muliplicity as well as the lower limit on W V i S reduce 
the background from 7*7* — > t + t~ events. The upper limit on W V i S cuts down the 
background from the e + e~ — > qq (7), e + e~ — > t + t~ and four-fermion processes. The 
comparison of the W vis distributions (Fig. 1) for the data and the Monte Carlo (MC) 
generated samples of events, described below, illustrates the effects of the W vis cuts. 

About 910k events are selected after application of the above selection criteria. 

3 Data Analysis and Results 

Monte Carlo samples of the various final states present in the data were generated for 
comparison with these data. The simulation of the process 7*7* — > hadrons was based 
on PYTHIA 6.143 [6] in which the description of the hadron production encompasses 
the processes described by the Quark Parton Model (QPM) (direct process), the Vec- 
tor Dominance Model (VDM) and the hard scattering of the hadronic constituents of 
quasi- real photons (resolved photon process). The MC sample of events used is 2.7 times 
larger than the data. The main background coming from the inclusive e + e~ — > qq (7) 
channel has been estimated from a PYTHIA 6.125 sample. The simulations of the 
e + e~ — > four-fermion, the 7*7* — > t + t~ and of the e + e~ — > t + t~ backgrounds were 
based on the EXCALIBUR [7], BDKRC [8] and KORALZ 4.2 [9] generators, respec- 
tively. The Monte Carlo generated events were then passed through the standard DEL- 
PHI detector simulation and reconstruction programs [5]. The same cuts were applied 
on the reconstructed MC events as on the data. 

The dN/dpT distribution of the charged particles of the selected events is presented 
in Fig. 2, for tracks with pseudo-rapidity \r)\ < 1 (77 = — In tan(0/2)) 1 , i.e. well 
measured tracks including TPC information. The expected Monte Carlo generated con- 
tributions, normalized to the data integrated luminosity are also shown. The data are 
well reproduced by the sum of the simulated samples of events for px > 1-6 GeV/c and 
the e + e~ — > qq (7) channel is the main contributor for pt > 12 GeV/c. There is a lack 
of data at pr < 1.6 GeV/c, becoming substantial at px < 1 GeV/c. This is caused by 
the trigger efficiency which was not accounted for in the Monte Carlo simulation and 
which is low for low pt tracks and low multiplicities [10]. For this reason, the following 
comparison with theoretical predictions is presented for p T > 1.6 GeV/c only. 

The differential da/dp T cross-section distribution of the inclusive production of charged 
hadrons in the process 7*7* — > hadrons has been obtained by subtracting the background 
contributions from the experimental dN/dpx data. The resulting distribution has been 

The angular selection of tracks (Table 1 and Figs. 2-5) is expressed in terms of |r)| cuts for comparison with published 
results [2,3]. 
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corrected, bin-by-bin, by a factor which is the inverse of the ratio of the numbers of 
reconstructed to generated tracks of 7*7* — > hadrons in Monte Carlo events. This ratio 
is of the order of 50-60% for 1.6 GeV/c < p T < 4 GeV/c and drops to about 20% for 
Pt > 10 GeV/c, the upper bound on W v i s being mainly responsible for the drop in 
efficiency on large pr tracks. The da/dpT distribution is shown in Fig. 3 for different sets 
of selection criteria as described below. The PYTHIA prediction is also shown. It agrees 
very well with the data for p T > 1.6 GeV/c up to large pt values. 

To study the systematic uncertainty coming from the selection criteria, we have varied 
them, in particular the W V i S upper limit and the track polar angle (9) cuts. A smaller 
upper bound of W V i S has the advantage of minimizing the background contributions espe- 
cially the e + e~ — > qq (7) one. Tracks at low polar angle are missing TPC measurements 
and are thus less well measured. On the other hand most contributing processes corre- 
spond to the emission of tracks peaked in the forward (backward) regions, in particular 
the e + e~ — > qq (7) and even more the 7*7* — > hadrons channels. Hence a tight (9) cut 
can reduce significantly the number of selected charged-particle tracks (N ch ) of a given 
event and consequently its computed visible energy W vis . Fig. 3 shows the da/dpr dis- 
tributions, calculated using tracks with \r)\ <1.5, for four sets of selection criteria varying 
the polar angle selection imposed on charged tracks and the cut on the visible invariant 
mass W V i S : 

1. 10° < 9 < 170° (\n\ < 2.4), 5 GeV/c 2 < W vis <20 GeV/c 2 ; 

2. 25° < 9 < 155° (\n\ < 1.5), 5 GeV/c 2 < W vis <20 GeV/c 2 ; 

3. 10° < 9 < 170° (\n\ < 2.4), 5 GeV/c 2 < W vis <35 GeV/c 2 ; 

4. 25° < 9 < 155° (\n\ < 1.5), 5 GeV/c 2 < W vis <35 GeV/c 2 . 

The spread of the measurements is relatively small for p T < 7-8 GeV/c but increases 
for high p T values where the e + e~ — > qq (7) dominates. The corresponding systematic 
uncertainty has been estimated as half of the spread of the four sets of measurements. 

The other source of uncertainty comes from the Monte Carlo modelling. It 
has been estimated by comparing the PYTHIA and TWOGAM [11] predictions for 
the 7*7* — > hadrons processes and PYTHIA and HERWIG [12] predictions for the 
e + e~ — > qq (7) process. It was found that the relative difference on the efficiencies cal- 
culated from the various generators depends on p? but never exceeds 10%. The corre- 
sponding uncertainty has been defined as half of the difference between two generator 
contributions. All systematic uncertainties have been added quadratically in Table 1. 

Table 1 gives the values of da/dp T as a function of p T) for the selection criteria de- 
scribed in section 2, the pseudo-rapidity ranges \n\ < 1 and \n\ < 1.5 and for pr > 1.6 
GeV/c where the event trigger efficiency is close to 100%. The first error is statistical and 
the second one is the overall systematic uncertainty. Fig. 4 shows the comparison of the 
da/dpr distribution for \rj\ < 1.5 with the NLO QCD prediction of [13]. The theoretical 
computation tends to be slightly lower than the measurements at high p T values although 
staying compatible with them within errors. 



4 Discussion of results 

Our measurement of the da/dpx cross-section of the inclusive production of hadrons in 
7*7* interactions appears to agree well with both PYTHIA and NLO QCD predictions. 

The L3 experiment has performed a similar analysis [2] and has observed that the 
Pt spectrum of charged hadrons is slightly below the PYTHIA MC prediction while the 
derived da/dpr cross-section considerably exceeds the NLO QCD prediction at high pt 
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values. We have repeated our analysis, adopting a "L3-like" set of selection criteria 
which, compared to ours, corresponds to a less tight W vis cut (W vis < 78 GeV/c 2 in- 
stead of 35 GeV/c 2 ) and a higher threshold of the total number of particles including 
neutrals (5 instead of 4). The looser W v %s cut has the effect of increasing significantly 
the e + e~ — > qq (7) background (see Fig. 1) which now dominates at large pt values. The 
resulting dN/dp T spectrum of charged particles for the "L3-like" events is presented in 
Fig. 5 together with the contributing channels. One observes an excess of data over the 
PYTHIA MC prediction. This disagreement between MC and data is likely to come from 
charged particles of background channels as these are introduced in much larger quanti- 
ties than charged particles from 7*7* — > hadrons, when the W V i S cut is relaxed up to 78 
GeV/c 2 , as can be checked by comparing Fig. 5 with Fig. 2. It legitimates, a posteriori, 
our < 35 GeV/c 2 cut to minimize the contamination of charged particles from 

background channels. 

The OPAL experiment has measured the differential da/dpx cross-section of the inclu- 
sive production of charged hadrons [3] for different intervals of W, the hadronic invariant 
mass corrected for detector effects. In the (10 GeV/c 2 < W < 30 GeV/c 2 ) interval, the 
cross-section is compatible with the NLO prediction. 

5 Conclusions 

The study of the inclusive charged hadron production in two-photon collisions has 
been carried out at the DELPHI detector at LEP II. Measurements of the inclusive 
single-particle pt distribution and of the differential inclusive da/dpx cross-section have 
been extracted. The differential inclusive da/dpr cross-section is found to be compatible, 
within errors, with the PYTHIA and NLO QCD predictions up to high p?, although 
systematic uncertainties limit the accuracy of the comparison in this region. It is shown 
that if cuts such as those used in [2] are applied, qq background dominates at large pr, 
making it difficult to draw conclusions on two-photon processes. 
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Table 1: Differential inclusive da/dpr of charged particles produced in 7*7* — > hadrons 
collisions, for \q\ <1, \r}\ <1.5 and pt > 1.6 GeV/c. The first error is statistical, the 
second is the systematic uncertainty. The data are background subtracted and corrected 
for detector inefficiency and selection cuts. 
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W vi5 , GeV/c 2 



Figure 1: W V i S distributions for the data and for the simulated 7*7* — > hadrons (medium 
cross-hatching), 7*7* — > r + r _ (second largest cross-hatching) , e + e~ — > qq (7),t + t~ 
(small cross-hatching) and e + e~ — > W + W~ (largest cross-hatching) events at a/s^T = 200 
GeV. 
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Pt, GeV/c 

Figure 2: p T distribution of charged particles of the selected sample of events, for 
1 77 1 <1 together with the Monte Carlo generated contributing processes: 7*7* — > hadrons 
(largest cross-hatching), e + e~ — > qq (7) (medium cross-hatching), e + e~ — > W + W~, t + t~, 
7*7* — > t + t~ (small cross-hatching). 
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Figure 3: Differential inclusive da/dpT distributions of charged particles with \r]\ <1.5, 
produced in 7*7* collisions, for different sets of initial selection criteria. (The lower limit 
of W V is was W V is > 5 GeV/c 2 ). The data are background subtracted and corrected for 
detector inefficiency and selection cuts. The line is the corresponding PYTHIA prediction 
for 7*7* — > hadrons. 
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Figure 4: Differential inclusive da/dpr distribution of charged particles produced in 7*7* 
collisions for |^| <1.5 and 5 GeV/c 2 < W vis < 35 GeV/c 2 . The original data sample used 
to extract this cross section included tracks with 10° < 9 < 170° < 2.4). The 
data are shown as points with statistical + systematical error bars. They are background 
subtracted and corrected for detector inefficiency and selection cuts. The line is the NLO 
QCD prediction of [13] for 7*7* — > hadrons. 
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Figure 5: pt distribution of charged particles of the event sample after application of the 
"L3-like" selection criteria, for \r]\ <1 and 5 GeV/c 2 < W vis < 78 GeV/c 2 , together with 
the Monte Carlo generated contributing processes: 7*7* — > hadrons (largest cross-hatch- 
ing), e + e~ — > qq (7) (medium cross-hatching), e + e~ — > W + W~, r + r _ , 7*7* — > r + r _ 
(small cross-hatching) . 



